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ABSTRACT 

We employ high-resolution dissipationless simulations of the concordance ACDM cosmology (Up = 
1 — Da = 0.3, h = 0.7, as = 0.9) to model the observed luminosity dependence and evolution of 
galaxy clustering through most of the age of the universe, from z ~ 5 to z ~ 0. We use a simple, non- 
parametric model which monotonically relates galaxy luminosities to the maximum circular velocity of 
dark matter halos (Umax ) by preserving the observed galaxy luminosity function in order to match the 
halos in simulations with observed galaxies. The novel feature of the model is the use of the maximum 
circular velocity at the time of accretion, V)))™ j for subhalos, the halos located within virial regions of 
larger halos. We argue that for subhalos in dissipationless simulations, U,((“ reflects the luminosity 
and stellar mass of the associated galaxies better than the circular velocity at the epoch of observation, 
U,(J°x . The simulations and our model L—Unax relation predict the shape, amplitude, and luminosity 
dependence of the two-point correlation function in excellent agreement with the observed galaxy 
clustering in the SDSS data at z ^ 0 and in the DEEP2 samples at z ~ 1 over the entire probed range 
of projected separations, 0.1 < rpj{h~^ Mpc) < 10.0. In particular, the small-scale upturn of the 
correlation function from the power-law form in the SDSS and DEEP2 luminosity-selected samples is 
reproduced very well. At z ^ 3 — 5, our predictions also match the observed shape and amplitude of 
the angular two-point correlation function of Lyman-break galaxies (LBGs) on both large and small 
scales, including the small-scale upturn. This suggests that, like galaxies in lower redshift samples, 
the LBGs are fair tracers of the overall halo population and that their luminosity is tightly correlated 
with the circular velocity (and hence mass) of their dark matter halos. 

Subject headings: cosmology: theory — dark matter — galaxies: halos — galaxies: evolution — 
galaxies:clustering — large-scale structure of the universe 


1. INTRODUCTION 

A generic prediction of high-resolution simulations of 
hierarchical Cold Dark Matter (CDM) models is that 
virialized regions of halos are not smooth, but con¬ 
tain subhalos — the bound, self-gravitating dark mat¬ 
ter clumps orbiting in the potential of their host halo. 
The subhalos are the descendants of halos accreted by a 
given system throughout its evolution, which retain their 
identity in the face of disruption processes such as tidal 
heating and dynamical friction. Their presence is in it¬ 
self a vivid manifestation of the hierarchical build-up of 
halo mass. 

In the CDM scenario, luminous galaxies form via cool¬ 
ing and condensation of baryons in the centers of the po¬ 
tential wells of dyk rnatter halos JWhite fc Reesl ll978t 
iFall fc Efstathioulll98nl: iBlumenthal et al.lll984fl . In the 
context of galaxy formation, there is little conceptual 
difference between halos and subhalos, because the lat¬ 
ter have also been genuine halos and sites of galaxy for¬ 
mation in the past, before their accretion onto a larger 
halo. We thus expect that each subhalo of sufficiently 
large mass should host a luminous galaxy and this is 
indeed sup ported by self-consistent cosmological simula¬ 
tions le.g.. iNagai fc Kravtsovll2?)?)^ . The observational 
counterparts of subhalos are then galaxies in clusters and 
groups or the satellites around individual galaxies. In 
this sense, we will use the term halos to refer to both 
distinct halos (i.e., halos not located within the virial 
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radius of a larger system) and subhalos. 

Although this general picture is definitely reason¬ 
able, it is not clear just how direct the relation be¬ 
tween halos and galaxies is. One may argue, for 
example, that the subhalos can be disrupted much 
faster than the more tightly bound ste llar system they 
host, leaving behind “orphan” ga laxies llGao et alJl2004 
iDiemand. Moore, fc Stadel 1200^ . At the same time, 
properties of surviving subhalos, such as maximum cir¬ 
cular velocity and gravitationally bound mass, are sub¬ 
ject to strong dynamical evolution a s they orbit within 
the potential of th e ir host halo (e . g..lMoore et alJll99fil: 
Kjyjhn eUaljJTggTlljjj yits^ iKravtsovetM] 

^OOdTlKazantzidiserahl 120041 . This makes the rela¬ 
tion betwee n subhalo properties and galaxy luminosity 
ambiguous l|Nagai fc Kravtso'^l200flll . because the latter 
may be less affected by dynamical processes but may 
evolve due to aging of the stellar populations after ram 
pressure strips the existing gas and the accretion of new 
gas is suppressed. The key question that we address in 
this paper is whether there is a one-to-one correspon¬ 
dence between populations of halos in dissipationless cos¬ 
mological simulations and galaxies in the observable uni¬ 
verse. As a test, we use comparisons of the predicted 
clustering of halos with the available observational mea¬ 
surements of galaxy clustering from z ~ 5 to the present. 

During the last decade, large observational surveys of 
galaxies both at low and high redshifts have tremen¬ 
dously improved our knowledge of galaxy clustering, its 
evolution, and the relation between the galaxy and mat- 
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ter distributions. A coherent picture has emerged in 
which bright galaxies are strongly biased with respect 
to the matter distribution at high redshifts llSteidel et al.l 


1998HGiavalisco et alJl99i4lAdelberger et al 


Ouchi et alJl2004bl [2001 iLee et alJl2005t iHamana et all 


2005(1 ■ and in which the bias decreases with time in such 
a way that the ampli tude of galaxy clus tering is only 
weakly evolving (e.g., lOuchi et alJl2004bl). as expected 
in hierarchical structu re formation llColm et alJ 1199^ 
iKaiiffmann et al.ll999ll . The bias is also in general scale- 
and luminosity-dependent. Bright (red) galaxies are 
more strongly cluster ed than faint (blue) galaxies both 
in the local universe (Norberg et al.l 121)02 iZehavi et al.1 
l2?M l2?M and references t herein) and in the distant 
past llCoil et a,1Jl20Ml200^ . 

A recent development is the detection of a departure 
from a pure power law in the two-point correlation func - 
tion of galaxies at z ~ 0 llZehavi et al.l 12004 1200511 . 
This departure is expected in general because the two- 
point correlation function is a sum of two separate con¬ 
tributions: the one-halo term, which arises from pairs 
of galaxies within a distinct dark matter halo, and the 
two-halo term, which arises fro m pairs of galaxies from 
two different distinct halos ('e.g.. lCoorav fc Shethll2003l . 
The one-halo contribution dominates on small scales, 
while at scales larger than the size of the largest viri- 
alized regions clustering is due to the two-halo term. 
The two terms are not generically expected to com¬ 
bine so as to give a power-law correlation function. 
The deviation of the correlation function from a power- 
law w as predicted t o be even stronger at higher red- 
shifts llZhengl 12004 iKravtsov et al.l T2 004bl). and this 
^has^ow been convincingjly confirined l| A(|(j^]l^(j^rg(j^r tj^t 

2005^ using gala xy samples identified wit h the Lyman- 


break technique 
Despite impressive advances in the amount and qual¬ 
ity of data on the galaxy distribution over a wide range 
of redshifts, the exact relation between dark matter ha- 


still rather uncertain fe.g.. 

Mo et al.lll999t IKolatt et all 

IT^ iWechsler et al.l Hfim 
tempts to connect the dar 
ies employ semi-analvtic m( 

). The most popula 
c halos and luminous j 
Ddeling fe.ET.. IWhite Ki. 

r at- 
ialax- 

Yen Id 

19911 iKauffmann et al. 

11994 lAvila-Reese et alJ 

199a 

Somerville & PrimackI 

19991 IGoleetalJ 11994 

200(1 

Groton et al.ll200,5l iBower et alJl2005ll. which uses 

phe- 


nomenological recipes for specifying when, where, and 
how galaxies form within dark matter halos, often in 
conjunction with high-resolution dissipationless simula¬ 
tions. A nother popular approac h is to use the halo model 
(see e.g. iCoorav fc Shetbll2002l and references therein), 
which, in its simplest form, specifies the probability dis¬ 
tribution for a particular halo of mass M to host a given 
number of galaxies N with specified properties, such as 
luminosity, color, etc. More complex halo models include 
the conditional luminosity functi on (CLF) approach 

which specifies the 


IjYa.ug et a,lJl2fln,4 ICoora,vll2r 
luminosity function for halos of mass M, and models 
which connect other features (e.g..lVale fc Ostrikerl2004 


iNevrinck. Hamilton, fc Gnedinll2005ll . 

While these approaches manage to capture the gen¬ 
eral observational trends, they usually employ a large 
number of free parameters, making it difficult to glean 


the relevant underlying mechanism(s) responsible for the 
agreement between the model and data. Most imple¬ 
mentations also do not provide a direct constraint on the 
relation between luminous components of galaxies and 
dark matter halos. In fact, assumptions are made about 
this relation in most cases. 

In this study we use galaxy clustering to address two 
straightforward issues. First, we question whether cur¬ 
rent data is consistent with a one-to-one correspondence 
between luminous galaxies and galaxy-sized dark mat¬ 
ter halos in cosmological simulations. Second, we ask 
whether the observed clustering is consistent with a sim¬ 
ple relation between galaxy luminosity and some prop¬ 
erty of its host halo. In the context of the second 
question, we investigate which halo properties are most 
closely related to galaxy luminosity. We show with the 
currently available data the answer to both of these ques¬ 
tions is yes, as we can reproduce luminosity-dependent 
clustering measurements at different redshifts from z ~ 5 
to z ~ 0 with a simple, non-parametric model relating 
galaxy luminosity to the halo circular velocity. 

Although a number of studies during the last 
decade have shown that galaxy clustering can 
be approximately matched by the clustering of 

dark matter halos in dissipationle ss _ si mulations 

(iGa.rlber'^ If 9911 iRraiuerd V Vilhimsenl 1 992L If 994hlal 
a)lm yrlberg fc CouchmanI 11994 lwSsle7et alJ 

Jc 


ilI98(_jGolfimd^Mj_Tl99ij^ 3Kraytsqv^/KK^ Il999t 

Kravtsov et al.ll2004bl: INevrinck et al.l l200'y the size of 
observational samples did not allow th orough tests of 
the gal axy-h alo relation. Fo r exam ple, IKravtsov et alJ 
<l2004bll and INevrinck et al.l ll2004l compare the two- 
point correlation function of halos and bright galaxies 
from the SDSS and PSGz surveys, respectively, and find 
good agreement on scales 0.1 < rj{h~^ Mpc) < 10. 
However, these studies do not attempt to match the 
clustering of fainter galaxies. iTasitsiomi et abl l|2004i 
assign luminosities to halos and compute the galaxy- 
mass correlation function, finding good agreement for 
two broad luminosity bins, after a reasonable amount of 
scatter was introduced into the relation. 

The current work extends these analyses by compar¬ 
ing results of very high-resolution simulations to the most 
current measurements of the two-point correlation func¬ 
tions over a wide range of luminosities and redshifts. The 
large size and wide luminosity range of observational 
samples allows us to test the relation between galaxy 
luminosities and properties of their host halos with un¬ 
precedented power. The novel feature of the simulation 
analysis we present is that for each halo and subhalo we 
track the evolution of its properties, such as mass and 
maximum circular velocities. As we show below, this is 
a key piece of information for reasons that are easy to 
understand. For distinct halos, the current circular ve¬ 
locity is a measure of their potential well built-up during 
evolution, and can therefore be expected to be tightly 
correlated with the stellar mass (or more generally the 
baryonic mass) of the galaxy the halo hosts. The cir¬ 
cular velocity of subhalos in dissipationless simulations, 
on the other hand, is a product of both mass buildup 
during the period when the halo evolved in isolation 
and tidal mass loss after the halo starts to orbit within 
the virialized regio n of a larger object and experience 
strong tidal forces l|Havashi et al.ll2fi0^ IKravtsov et alJ 
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l2004at iKazantzidis et alJ 120041) . The stellar component 
of the galaxies, which should be more tightly bound than 
halo dark matter, should not be significantly affected by 
tidal forces and can stabilize the mass distribution (and 
hence Umax) in the inner regions. We can therefore ex¬ 
pect that luminosity and stellar mass of galaxies hosted 
by halos in dissipationless simulations should be corre¬ 
lated with the subhalo mass or circular velocity, , 
at the epoch of accretion, rather than with its current 
value. This is borne out by cosmological simulations, 
which include gas dynam ics, cooling, and star formation 
(iNag-ai fc Kravtso'^l200fijl who show that selection using 
Knax results in subhalo properties similar to the selection 
based on stellar mass of galaxies subhalos host. One can 
therefore argue that a reasonable approach is to relate 
galaxy luminosity to the current halo circular velocity for 
distinct halos and to the circular velocity at accretion for 
subhalos. The main result of this study is that this simple 
model reproduces the luminosity-dependence of galaxy 
clustering at different epochs with rem arkable, and per¬ 
haps s urprising, accuracy. We note that iVale fc Ostriked 
11200511 have recently presented and used a semi-analytic 
model for subhalos, which employs a similar approach 
to luminosity assignment, except that they use the total 
bound halo mass instead of circular velocity. 

The paper is organized as follows. In we briefly de¬ 
scribe the simulations, halo finding algorithm, and the 
method for tracking the evolution of halos. details 
our method for relating halos to galaxies, and elabo¬ 
rates upon our motivation for using the basis 

for the luminosity assignment for subhalos. In m we 
compare observational clustering results over the red- 
shift interval 0 < z < 5 to the clustering of halos in 
dissipationless AGDM simulations. The halo occupa¬ 
tion distribution implied by this model is described in 
33 In 0 we discuss the implications of our results and 
in 33 we summarize our main conclusions. Through¬ 
out this paper we assume a AGDM cosmology with 
{nm,nA,h,as) = (0.3,0.7,0.7,0.9). 

2. THE SIMULATIONS 

The simulations used here were run usi ng the Adaptive 
Refinement Tree (AR T) Wbody code llKravtsov et al.l 
Il997t IKravtso-^IlQQOU . The ART code implements suc¬ 
cessive refinements in both the spatial grid and temporal 
step in high density environments. These simulations 
were run in the concordance flat AGDM cosmology with 
ftm = 0.3 = 1 — Ua, h = 0.7, where ftm and Da are the 
present-day matter and vacuum densities in units of the 
critical density, and h is the Hubble parameter in units of 
100 km s“^ Mpc“^. The power spectra used to generate 
the initial conditions for the simulations were determined 
from a direct Boltzmann code calculation (courtesy of 
Wayne Hu). We use a power spectrum normalization of 
as = 0.90, where as is the rms fluctuation in spheres of 
8h~^ Mpc comoving radius. 

To study the clustering properties of dark matter halos 
over a range of scales, we consider two simulations of the 
above cosmology The first simulation, L80, followed the 
evolution of 512^ particles in a 80/i“^ Mpc box. The L80 
simulation has a particle mass of nip = 3.16 x 10® h~^MQ 
and peak force resolution /ipeak = 1-2 h~^ kpc. This 
is the simulation to which we make most comparisons 
with observations. The second simulation we consider is 



100 1000 
v:;,“ (kms-‘) 


Fig. 1.— Top panel: Cumulative velocity function for all halos 
identified in the L80 simulation at various redshifts, in units of 
Mpc~^. Bottom panel: Fraction of subhalos as a function of 
redshift and maximum circular velocity at the time of accretion, 
KnS • We truncate the curves where A^sub < 10 because in that 
regime poisson noise washes away any useful information. The 
arrow delimits our nominal completeness limit. 


denoted L120 and was run with 512® particles in a 120 
Mpc box, resulting in a particle mass of nip = 1.07 x 
10® h~^MQ and peak force resolution /ipeak = 1-8 h~^ 
kpc. This simulation thus has a larger particle mass and 
somewhat lower spatial resolution compared to the L80 
run. We use this simulation to obtain better statistics for 
the correlation function of rare (i.e., massive) objects. 

2.1. Halo Identification, Classification, and 
Construction of Merger Trees 

Our analysis requires detailed dynamical knowledge of 
not only distinct halos, i.e. halos with centers that do not 
lie within any larger virialized system, but also subhalos 
which are located within the virial radii of larger systems. 
Note that the term “halo” (e.g., the halo occupation dis¬ 
tribution) usually refers to what we call distinct halos in 
this work. 

We identify distinct halos and the subhalos within 
them using a variant of the Bound Density Ma xima 
(BDM) halo finding algorithm lIKlvoin et al.lll99^ . De- 
tai ls of the algorithm and p arameters used can be found 
in iKravtsov et al.l ll2004 we briefly summarize the 
main steps here. All particles are assigned a density us¬ 
ing the smooth algorithm^ which uses a symmetric SPH 
smoothing kernel on the 32 nearest neighbors. Starting 
with the highest overdensity particle, we surround each 
potential center by a sphere of radius rgnd = 50/i“® kpc 
and exclude all particles within this sphere from further 
search. Hence no two halos can be separated by less than 

^ To calculate the density we use the publicly available 
code smooth: http://www-hpcc.astro.washington.edu/tools/ 
tools.html 
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?'find- We then construct density, circular velocity, and 
velocity dispersion profiles around each cent er, iteratively 
remov ing unbound particles as described in lKlvoin et al.l 
(TT^ . Once unbound particles have been removed, we 
measure quantities such as Knax = \/GM[< r)/rl^ax) 
the maximum circular velocity of the halo. For each dis¬ 
tinct halo we calculate the virial radius, defined as the 
radius enclosing overdensity of 180 with respect to the 
mean density of the Universe at the epoch of the out¬ 
put. We use this virial radius to classify objects into 
distinct halos and subhalos. The halo catalog is com¬ 
plete for halos with more than 50 particles which corre¬ 
sponds, for the L80 simulations, to a minimum halo mass 
of 1.6 X lO^O/i-i Mq. 

For subhalos we also tabulate l^ax j maximum cir¬ 
cular velocity at the time when a subhalo falls into a 
distinct halo. In order to tabulate this quantity we rely 
on merger trees generated for these simulations, which 
track the histories of both distinct halos and subhalos. 
A detail ed description o f the merger tree construction is 
given in lAllgoo'dr ^|2005^ . The merger trees we use here 
follow halo evolution through 48 timesteps between z = 1 
and the present for the L80 box and 89 timesteps for the 
L120 box. For each subhalo, we use the merger trees to 
step back in time until the subhalo is no longer identified 
as belonging to a larger system. We then define to 
be Vnax of the subhalo at that time. 

In the top panel of Figure ^ we show the cumulative 
velocity function for all identified halos from z = A to 0. 
This figure quantifies the relation between 1^“ 
we will use these two quantities interchangeably to define 
our halo samples throughout the paper. In the bottom 
panel we show the corresponding cumulative fraction of 
subhalos as a function of time. The figure shows that 
the subhalo fraction is a weak function of circular veloc¬ 
ity at all epochs. There is a weak trend for a smaller 
subhalo fraction among halos with larger lUJ™ ■ There 
is a stronger trend of increasing subhalo fraction with 
decreasing redshift. 


3. CONNECTING GALAXIES TO HALOS 

In this section we motivate and describe our model for 
associating galaxies with dark matter halos. We make 
this connection by relating galaxy luminosity and a phys¬ 
ical property of dark matter halos, for which we choose 
KSax foi” distinct halos and 1^“ for subhalos. As dis¬ 
cussed above, for subhalos denotes the maximum 
circular velocity at the time the subhalo was accreted. 
The maximum of the circular velocity profile, Umax) is 
a measure of the depth of the halo potential well and 
is expected to correlate strongly with stellar mass of 
the galaxies, as implied by the Tully-Fisher and Faber- 
Jackson relations. At the same time, the definition of 
Umax in simulations is unambiguous both for distinct ha¬ 
los and subhalos, which is not the case for the total mass, 
as different operational definitions are used by different 
authors. It should be noted that Umax used here will not 
correspond directly to Umax observed in, for example, the 
Tully-Fisher relation because dissipationless simulations 
do not take into a ccount the effect of bar yon condensa¬ 
tion on Umax l'e.g.. lBhimenthal et ahlllQ^ . 

The use of subhalos is the novel feature 



O.I 1.0 0.1 1.0 10.0 

fp (h-' Mpc) 


Fig. 2.— Projected two-point correlation function at 2 : 0 

comparing the effects of selecting on (solid lines) versus 

Knax (dashed lines) at four different number density thresholds 
(labeled in the top right corner, in units of Mpc“^). While 
there is a slight increase in the correlation function on large scales 
when using rather than VjJJSx ? the difference is much stronger 
on small scales. The difference between and is due to 

max max 

the tidal stripping of subhalos which have fallen into larger sys¬ 
tems, hence correlation functions will be most strongly effected on 
small scales. 


of our modeU. As we discussed in the introduc¬ 
tion, the motivation for this is fairly straightfor¬ 
ward. While Vmax decreases due to tidal stripping as 
a halo falls through a larger halo (jHavashi et a,1.ll2flfl3l 
iKravtsov et al.ll2004all. one can expect that the stellar 
component of galaxies within these halos will not be af¬ 
fected appreciably since stars are concentrated near the 
bottom of t he halo potential well an d are more tightly 
bound le.g.. iNagai fc KravtsovI 1^00511 . Hence we argue 
that, for galaxies associated with subhalos, luminosity 
should correlate more strongly with than the U,Sax 
affected by dynamical evolution. Therefore, throughout 
the rest of the paper, unless explicitly stated otherwise, 
the maximum circular velocity, Umax , will be assumed to 
mean: 

,r f T,Saxj subhalos 
max = I distinct halos 

Although it is not clear how accurate this assumption is 
in detail, we show that it provides a considerably bet¬ 
ter match to observed galaxy clustering compared to the 
uniform selection using U^jJ™ for both subhalos and dis¬ 
tinct halos. Note also that the use of circular velocities 
before accretion can also help explain t he abundance of 
faint dwarf galaxies in the Local Group llKravtsov et alJ 

mi. 

In order to assign luminosities, we assume a mono¬ 
tonic relation between galaxy luminosity and Umax and 

^ A s we prepared this paper for publication, |YaJe_^,22llik£3 
ll2()n^ submitted a paper in which they employ a semi-analytic 
model for subhalos and a similar approach to luminosity assign¬ 
ment, except that they use the total bound halo mass instead of 
circular velocity. 
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require that the L—Vmax relation preserves the galaxy 
luminosity function (LF). Specifically, we use the follow¬ 
ing equation: 


Li) — 7Z/j(> (2) 

where Ug and Uh are the number density of galaxies and 
halos, respectively. We stress again, that Umax in the 
above expression is equal to U,^ax subhalos and to 
Knax foi' distinct halos (for which the accretion epoch 
is undefined). For each Li we find the corresponding 
Umax.i such that the above relation is satisfied. The main 
assumption is therefore that there is a monotonic rela¬ 
tion between galaxy luminosity and Umax ■ The model 
makes no further assumptions and is completely non- 
parametric. This should be kept in mind when we com¬ 
pare predictions of this model to observed galaxy clus¬ 
tering. Note that we do not take any possible scatter in 
the L— Unax relation into account in this model. 

With the L—Unax relation in hand, comparing obser¬ 
vational clustering statistics to the model predictions is 
straightforward: we simply compute the desired statistic 
for the halos with assigned luminosities corresponding to 
the observed sample luminosity range. This method does 
not currently treat other galaxy properties such as color, 
although it could conceivably be extended to include such 
properties. We have not included the possibility of “or¬ 
phaned galaxies”, i.e. galaxies without any associated 
subhalos. We discuss the issue of orphans in detail in m 

Figure 121 shows the effect that selection of halos using 
U,0“x rather than for subhalos has on the projected 
two-point correlation function iOp for different number 
densities, at z = 0. As expected, the effect is most sig¬ 
nificant on small scales, where the subhalo contribution 
dominates, and the difference between U^™ and U0°x is 
greatest. 

In Figure O we show the effect of selecting halos 
according to U,0“x and U,))™ on the galaxy-mass cross¬ 
correlation function, ^gm- The small bump in the sam¬ 
ple selected using UO™ in the top panel is due to the fact 
that U^“ samples in general have more satellites than 
Knax samples, and the satellite contribution to ^gm peaks 
near Vp = 0.5 h~^ Mpc. The “bump” is smaller in the 
lower panel, where the Umax threshold is much higher, be¬ 
cause the number of satellites is nearly the same between 
the Umax s-iih T)0ax Selected samples in this case. 

At higher redshifts the differences between samples 
selected using U^™ T,nax small.Figure 0 

shows the projected correlation function for the USax” 
and Umax -selected samples at four redshifts. The sam¬ 
ples are constructed to have a fixed number density 
n = 1.5 X Mpc^. The situation is similar for 

other number densities. Already by 2 : ^ 1 the effect of 
selection is quite small. The difference on small scales 
for the z ~ 3 and z ^ 4 samples is not statistically sig¬ 
nificant.^ 

Note that the designation “now” refers to the time of observa¬ 
tion, not to z = 0. 

^ Note that at high 2 : the correlation function at the smallest 
scales is somewhat higher for the -selected samples, which 

seems counter-intuitive. We believe that this is a small artefact 
of the halo finding algorithm. At higher redshifts, the halos are 
smaller and subhalos are typically at smaller distances from the 
host center. At small radii the removal of unbound particles is 
more difficult as the halos are located near the bottom of the po- 



0.1 1.0 10.0 
Tp (h-‘ Mpc) 


Fig. 3.— Comparison of the galaxy-mass cross-correlation func¬ 
tion for halos selected using (solid lines) and (dashed 

lines) at two different circular velocity thresholds (labeled in the 
top right corner of each panel, in units of km s~^). 



0.1 1.0 0.1 1.0 10.0 
rp(h-‘ Mpc) 


Fig. 4. — Comparison of the projected two-point correlation func¬ 
tion for halos selected using VmSc (solid lines) and (dashed 

lines) at four different redshifts, for a fixed number density, n = 
1.5 X Mpc“^. This figure clearly shows that, while using 

KnS over results in a large difference at low redshift, it has 

very little impact at higher redshifts. The trend is similar for a 
wide range of number densities. 


tential well. The value of such cases can be overestimated 

producing a boost in the number of subhalos above a given thresh¬ 
old value, and boosting the correlation function somewhat. Note, 
however, that the effect is small and the difference between corre¬ 
lation functions for VJ^ax " bJJax -selected samples is less than 
2cr. 
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We believe that this effect can be attributed to the 
distribution of accretion times for subhalos at each red- 
shift: at low redshift, subhalos have a wide distribution 
of accretion times, and hence a large number of subha¬ 
los have had time to experience significant tidal strip¬ 
ping, while at higher redshifts the distribution of accre¬ 
tion times rises sharply near the epoch of observations. 
This is because both the accretion and disruption rates 
are high at high redshifts. The accreted halos do not sur¬ 
vive for a prolonged period of time, so that at each high-z 
epoch most of identified subhalos are recently accreted 
objects, which are yet to experience significant tidal mass 
loss. 

Since we are only computing for subhalos which 
have survived to the current epoch, one may worry that 
we are neglecting a significant population of subhalos 
with a sizeable E^ax that are not present in the halo cat¬ 
alog at the current epoch. Of such a population there 
can only be two fates: either the object was at some 
point physically disrupted or it has simply fallen below 
the resolution limit. We have used the merger trees to 
find all subhalos which have ever fallen into a distinct 
halo and have tabulated their and the ratio in mass 
between the subhalo and distinct halo, at the time of 
accretion. For a wide range in thresholds, the dis¬ 
tribution of mass ratios is strongly peaked between 0.1 
and 1.0. This implies that dynamical friction has caused 
the subhalo to merge with the distinct halo on the order 
of a dynamical time, and suggests that the majority of 
these subhalos have in fact physically merged and should 
not have survived. The absence of such missing subhalos 
also implies that in our simulations there should be no 
“orphan” galaxies lIGao et al ] |20M see § ini for further 
discussion of this issue). 

4. GALAXY CLUSTERING FROM Z ~ 5 TO THE PRESENT 

In this section we compare clustering statistics of halos 
to recent observations of the galaxy two-point correlation 
function over the redshift interval 0 < z < 5. The ob¬ 
served clustering statistics we compare to are ujp(rp), the 
projected two-point correlation function, and oj{0), the 
angular two-point correlation function. We estimate ujp 
by integrating the real space, three-dimensional correla¬ 
tion function, ^(r), computed in the simulation, along 
the line of sight: 


Wp(rp) — 2 / 

Jo 

where the comoving distance r has been decomposed into 
perpendicular (r^) and parallel to the line-of-sight (ry) 
components. In practice, the integration in Eqn. |3 is 
truncated at some finite scale: we truncate at 40/i“^ 
Mpc for SDSS galaxies f U4.1 ll and 2Qh~^ Mpc for DEEP2 
galaxies (S3, as is done in the data. Since the simula¬ 
tion box size is only 80h~^ Mpc, the measurement of ^(r) 
is not reliable for r > O.lLbox ~ 8h~^ Mpc. To extrap¬ 
olate ^(r) to larger scales, we use the two-point 

correlation function of dark matter multiplied by the 
linear bias of ^(r) measured over 4 < r/(h“^Mpc)< 8. 



® We derive the dark matter correlation function fromUimgower 
s£ecHum provided by the publicly available code of |Smitli^t_aJJ 
(2223), which is more accurate than the popular Peacock and 
Dodds prescription. 


Generating uj{9) from ^(r) without assuming ^(r) to be 
a power law is somewhat more involved. With knowledge 
of the redshift distribution, N{z), of the sample, oj{0) can 
be derived via the Limber transformation: 


U}{0) = 


/“ dzN^{z) dx^iy/[D^{z)0]^ + x^)/Rh{z) 
lIo^dzNizW 


( 4 ) 

where Dm{z) is the proper motion distance and Rh{z) 
is the Hubble radius at redshift z. As for Wp, the integral 
over ^(r) is in practice truncated at a finite scale; we 
integrate to 50 h~^ Mpc and note that the resulting oj{0) 
is not sensitive to this particular truncation scale. 


4.1. Clustering at z ~ 0 

The SDSS l|York et al.ll2000t lAbazaiian et al.ll2004D is 
a large phot ometric and spectros copic survey of the lo¬ 
cal Universe. iZehavi et al.l ll2005D have measured the lu¬ 
minosity and color dependence of Wp{rp) for ~ 200,000 
SDSS galaxies over « 2500 deg^ with z < 0.15. As men¬ 
tioned in m assigning Umax to galaxy luminosity while 
preserving the observed luminosity function (LF) results 
in a unique T—Umax relation. In order to make the as- 
si gnment we use t he SP SS lur ninosity function presented 
in iBlanton et al.l 1200.11) . with iSchechted 119761) parame¬ 
ters in the r-band M* — 51og/i = —20.5, a = —1.05, and 
(j)* = 1.5 X Mpc^. It is then straightforward to 

compare the observed luminosity dependence of both the 
small and large scale clustering of SDSS galaxies to our 
model. 

The results for luminosity threshold samples (L > Lth) 
are shown in the l eft pa nel of FigureEl where we compare 
the IZehavi et all 1200,51) results to the clustering of halos 
corresponding to the range of galaxy luminosities in each 
sample. For the three halo samples with n = 6 x 10~^h^ 
Mpc-3, n = 1.5 X Mpc-3 and n = 2.8 x IQ-^h^ 

Mpc“^ we use the L80 simulation, while for the halo 
sample with n = 1.1 x Mpc“^ we use the L120 

simulation in order to improve statistics. See Table^for 
details of the SDSS samples used here. The agreement is 
excellent over all scales. We find similar agreement when 
Wp{rp) is measured in differential, rather than integral, 
luminosity bins. It is critical to realize that the agree¬ 
ment on scales Vp ^ lh~^ Mpc is due to our luminosity 
assignment scheme using . The luminosity assigned 
using for subhalos would result in a significant 

under-prediction of amplitude of Wp at rp < lh~^ Mpc, 
especially for fainter samples (see FigureEj. 

The halo occupation distribution (HOD) specifies 
the distribution of the number of galaxies within 
a (distinct) halo of mass M, P(N\M). It has be- 


come a doi 

l.lins- et all 

Dular tool for interoreting galaxv clustering 
119981: ISeliakI 120001 IScoccimarro et all 

20011 

iBullock. Wechsler. Kr. Somervilld 

20021 

Ya,n. Ma.dg’ 

wick. Kr. Whitd 1200,11 iBerlind et al.1 

2(m 

Zehavi et al.11200,51: IKravtsov et all I2004bl: IZhend 

2004 

Abazaiian 

et all 1200,5l iTinker et al.1 l200,5l). which re- 


quires that the first two moments of P{N\M) be 
specified to calculate two-point clustering. In the right 
panel of Figure |31 we show the first moment of this 
distribution, the average number of galaxies within 
a (distinct) halo of mass M, {N{M)), for the four 
halo samples which correspond to the four luminosity 
threshold SDSS samples in the left panels. As expected. 
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Fig. 5.— Left: Comparison between the SDSS projected correlation function (points) and the correlation function derived from halos 
(solid lines) for various luminosity threshold samples. For comparison we include the correlation function of dark matter particles (dotted 
lines) at the median redshift of the sample. Right: The first moment of the halo occupation distribution (HOD) for the four halo samples. 
For all four samples, the gradual roll-off at small mass is due to scatter in the Hmax -mass relation. The fan (dotted lines) corresponds to 
slopes of 0.4, 0.7, and 1.0. 


the halo samples corresponding to brighter galaxy 
samples reside preferentially in more massive distinct 
halos. The halo sample corresponding to the brightest 
galaxies {Mr — 51og/i < —21) rarely has more than one 
halo per distinct halo. All three halo samples display a 
gradual roll-off in {N{M)) at low mass which is simply 
due to scatter in the Umax-mass relation, as we select 
samples using Umax , but plot as a function of mass. See 
for a more detailed discussion of the HOD associated 
with this model. 

The good agreement between the observed galaxy cor¬ 
relation function and samples of halos with our L—Umax 
model, over a range of luminosities and scales, suggests 
that the luminosity dependence of galaxy clustering is 
due primarily to how galaxies form within dark matter 
halos. This implies that galaxy properties vary as a func¬ 
tion of larger scale environment only insofar as the halos 
in which the galaxies reside vary. 

4.2. Clustering at z ~ 1 

Th e DEEP2 Galaxy Redshift Survey llDavis et al.l 
l2004^ has gathered optical spectra for ~ 50, 000 galaxies 
at z ~ 1 using the DEIMOS spectrograph on the Keck H 
10-m telescope. The survey, recently completed, spans a 
comoving volume of ~ Mpc^, covering 3 deg^ over 

four widely separated fields. We use the DE EP2 H-band 
luminosity function of lWillmer et al.l (120051) to compute 
the L—Vmax relation at z ~ 1. A Schechter fit to the 
overall luminosity function yields — 51og/i = —20.73 
and (j)* = 8.7 x Mpc^ with a fixed at a = —1.30. 

A detailed comparison has shown that these values are 
consistent with ot her estimates of th e global luminosity 
function at z ~ 1 llFaber et al.ll2005j) . 

The projected two-point correlation function, utp(rp), 
has been measured for DEEP2 galaxies as a function of 



0.1 1.0 0.1 1.0 10.0 
(h-‘ Mpc) 


Fig. 6.— Projected two-point correlation function at 2 ~ 1 for 
DEEP2 galaxies (solid circles) and halos (solid lines), at four differ¬ 
ent luminosity thresholds. We include jack-knife errors, computed 
using the eight octants of the simulation cube, on the model pre¬ 
diction for the brightest sample to demonstrate that they agree 
within Icr. The excellent agreement on all scales for these four 
samples suggests that luminosity-dependent clustering is a result 
of two effects: a simple relation between galaxy luminosities and 
dark matter halos, and the spatial clustering of the halos. For com¬ 
parison, we include the correlation function of dark matter particles 
(dotted lines). 


lumino sity and co lor lIGnil et a1JI2f)f)ll2nn5blal) . In ad¬ 
dition, iGoil et M I ll2nn5bD has estimated the two-point 
cross correlation between galaxies and groups, and be- 
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tween group cent ers and group gala xies, based on the 
group catalog of lOerke et al.l ll200fill . For the follow¬ 
ing comparisons we use the most recent me asurements 
of ujr,( r„) derived from the completed survey llCoil et al.l 

l2(M . 

Figure El compares the projected two-point correlation 
function for DEEP2 galaxies in four luminosity threshold 
samples to the clustering of corresponding dark matter 
halos. For the Mb — 51og/i < —19.0, —19.5, and —20.0 
samples we use the L80 simulation, while for the Mb — 
51og/i < —20.5 sample we use the L120 simulation to 
improve statistics. Slight discrepancies on small scales 
for the Mb — 51og/i < —20.5 sample may be attributed 
to cosmic variance and poisson noise in a sample of this 
number density, and in fact our smaller L80 box provides 
a slighly better match to the data. Overall the agreement 
is excellent on all scales for all four samples^. 

We would like to stress again that this remarkable 
agreement is achieved using the halo distribution in dis¬ 
sipationless simulations using a simple, non-parametric 
relation between galaxy luminosity and halo circular ve¬ 
locity. The luminosity-dependent bias at 2 ~ 1 hence 
seems to be driven entirely by the fact that brighter 
galaxies reside in more massive halos, with the corre¬ 
spondence between halo and luminosity determined by 
matching the observed luminosity function to the dark 
matter halo velocity function. 

4.3. Clustering at z > 2 

Very little was known about the clustering properties 
of galaxies at z > 2 until the advent of the Lyman- 
break technique I Steidel fc Hamiltonlll99,3l: ISteidel et al.1 
IT^ IT^ i2nn.C ~ This technique allows the identifi¬ 
cation of high-redshift Lyman-break galaxies (LBGs) by 
optical photometry alone, using information about the 
well-defined region in color-color space that these ob¬ 
jects occupy. While simple color-color cuts allow one 
to gather large numbers of LEG candidates with relative 
ease, it should be kept in mind that this technique is not 
perfect. Accurate completeness numbers are difficult to 
estimate, but it is believed that the Lyman-break tech¬ 
nique successfully identifies 80 — 90% of real LBGs (spe¬ 
cific completeness numbers depend on limiting apparent 
magnitude, dataset — e.g., ground observations or the 
HST, and sample definition, among other unknowns, and 
is often estim ated via Monte Carlo simulations of artifi¬ 
cial LBGs: seel^elberger et al.l20d^lOuchi et al.l2004at 
iLee et af]l2005lb 

The level of contamination, or the fraction of false posi¬ 
tives, is more important when considering the clustering 
of LBGs because such objects can artificially dilute or 
enhance the observed signal. At z ~ 3 the main source 
of contamination is Galactic stars (4%) and high-redshift 
active galactic nuclei f3 %b as determined by extensive 
spectroscopic follow-up llSteidel et alJl200^ . At higher 
redshifts the situation is less certain, as there has to 
date b een much less spectroscopic follow-up. lOnchi et al.1 
ll2005ll estimate the contamination for their z ~ 4 sam¬ 
ple based on Subaru data at 5% from spectroscopic 
follow-up of 63 LBG candidates. LBG candidates ex¬ 
tracted from the GOODS survey avoid Galactic star con¬ 
tamination thanks to high angular resolution of the HST 

^ Booyah! 



Mass (h'‘ M®) 


Fig. 7.— Comparison between z ^ LBGs and dark matter 
halos in simulations. Top panel, the projected correlation function 
measured from 700 spec troscopically confirmed LBGs (solid circles, 
lAdelherger et ?^l2nn.?ll . compared to halos at the same number 
density, n = 4 x Mpc“^ (solid line) and halos twice as 

numerous, n = 8 x 10~^h^ Mpc~^ (dashed line). Bottom panel. 
The number of halos per distinct halo for the two halo samples. 
The fan (dotted lines) corresponds to slopes of 0.4, 0.7, and 1.0. 


llLee et al ] |200l . although other sources of contamina¬ 
tion remain unquantified. 

In this section all quoted number densities have been 
completeness and contamination corrected. If the com¬ 
pleteness is well estimated, then using the corrected num¬ 
ber density should mitigate any incompleteness effects. 
However, even if the number density is contamination 
corrected, in order to fairly compare the clustering of 
LBGs to dark matter halos we must include the clus¬ 
tering properties of these contaminants. For simplicity, 
when estimating the effects of contamination we will as¬ 
sume the contaminants to have a random distribution 
on the sky. (We need not include additional clustering 
effects in the completeness case because in this case the 
missing objects are LBGs, which are assumed to have 
the same clustering properties as the identified LBGs, 
but note that this assumption would break down if the 
LBG completeness was a stro ng function o f lumin osity.) 

Since the early result llSteidel et al.1 1199811 that 
Lyman-Break Galaxies are a strongly clustered pop¬ 
ulation, there have been numerous attempts to 
use their clustering properties to connect these 
galaxies to their host dark matter halos (e.g.. 


Wechsler et al.ll998HJing fc Sutoll998HAdelberger et al. 

19981 iKatz et al.l 119991: iMo et al.1 

19991 1 Wechsler et al 

20011 IGia.valisco V DickinsonI 12001 

; IBullock et all 120021 


Scannapieco fc Thacked 1200^ . There have been two 


popular explanations for the properties of these galax¬ 
ies: some have speculated that they are a quiescent star- 
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Fig. 8.— Left: Two-point correlation function of 2 ~ 4 LBGs derived from Subaru data (solid circles. FOnchi et a] .1200^ compared to the 
correlation function for halos at the same number density (solid lines). The four panels correspond to four different apparent magnitude 
limits, and hence different number densities (see Table U for a summary of the Subaru data). The dotted line is uj{0) for dark matter 
particles at 2 ~ 4. Right: Average number of member galaxies for halos corresponding to the four Subaru samples at 2 ~ 4. The lines 
correspond to the Subaru samples, from left to right, i < 26.5, i < 26.0, i < 25.5, and i < 25.0. The fan (dotted lines) corresponds to slopes 
of 0.4, 0.7, and 1.0. 


forming populat i on of the most massive galaxies (e.g. 

Coles et al T, al JlQ^lGiavalisco fc Dickinsoril 

200 ID . while others have suggested that they are 
more relate d to temporal events such as merger-driven 
starbursts llLowenthal et alJ IIOOTI iKolatt et alJ Il999t 

iSomerville et alJ I2001I iScannapieco k. Thack^ I200,T 

The extent to which either of these scenarios can be ruled 
out by the clustering data has been a matter of some de¬ 
bate, depending on th e detailed assumption s that were 
made in each case (see iWechsler et al.ll2001l for a review 
of these issues). 

Our results for LBGs at z ^ 3 are shown in Figured 
In the top panel we plot the two-point projected correla- 
tion function for 700 sp ectroscopically confirmed LBGs 
ijAdelberger et al.ll200,‘lll with an estimated number den¬ 
sity of n = 4 X 10~^h^ Mpc“^ at z = 2.9. Note that this 
sample does not suffer contamination problems because 
these LBGs are spectroscopically confi rmed. We plot a;„ 
for ha los in units identical to those in lAdefberger et abl 
ll2?ifl!il where rmax is the line-of-sight distance through 
which we count projected pairs). We measure 0 Jp{rp) for 
halos at two different number densities to illustrate that 
a factor of two uncertainty in the LB G number density 
(a lar ger uncertainty than is quoted in lAdelberger et idl 
l2?in!ili does not alter these results. The figure shows that 
the agreement between the clustering of halos and LBGs 
is quite good. 

In the bottom panel of Figured we show the halo occu¬ 
pation of galaxies for the halo samples we associate with 
z ~ 3 LBGs. As expected, the halo sample which has a 
higher number density (dashed lines), and a correspond¬ 
ingly lower Umax threshold, has a lower cutoff in {N{M)). 
The measured halo occupation implies that most distinct 
halos are host to a single LBG. {N{M)) rises above 2 
only for the most massive halos (> which 



1 10 1 10 100 
0 (arcsec) 


Fig. 9.— Two-point corre lation function o f LBGs measured using 
GOODS data (solid circles. Fbee et af 11200^) . compared to the cor¬ 
relation function of halos at the same number density (solid lines). 
We also show the effect of a 10% contamination in the identifica¬ 
tion of LBGs (dashed lines). Top panels-. B-band dropouts which 
correspond to 2 4, at apparent magnitude limits of mz < 26.0 

(left) and tUz < 26.5 (right). Bottom panels: V-band dropouts 
corresponding to 2 ~ 5 at apparent magnitude limits, mz < 26.0 
(left), and mz < 26.5 (right). See Table [Tl for a summary of the 
GOODS data used here. 


are very rare at z ~ 3. 

At still higher redshifts the clustering of LBGs has 
recently been measured in two separate datasets: by 
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lOi^i et alJ 1 I 2 OO 5 D using Subaru data and bv iLee et al.1 
1120051) using the HST GOODS survey. These authors de¬ 
tect a strong departure from a power law in the angular 
two-point correlation function of LBGs. When parame¬ 
terized in the framework of the halo model, this excess 
small scale {9 < 10") power is attributed to multiple 
galaxies within a distinct halo, i.e. the ‘1-halo’ term. We 
now compare our L—l)nax model to these data. 

The left panel of Figure IH1 c ompares the a ngula r corre¬ 
lation function of LBGs from lOuchi et al.l i2005l) to the 
angular correlation function of halos for four different 
apparent magnitude thresholds, at 2 ~ 4. If the identi¬ 
fied LBGs occupy a relatively narrow range of redshifts, 
as expected, then these apparent magnitude cuts should 
correspond to absolute magnitude limits. The agreement 
is again remarkably good on all the scales we probe. Our 
non-parametric L—14iax model captures the luminosity 
dependence of LBG clustering, on both large and small 
scales, correctly predicts the scale of the small-scale up¬ 
turn in ui{0), and the fact that this scale decreases for 
increasing number density. The latter trend is due to 
a decrease in the typical mass (and size) of the distinct 
halos hosting LBGs as the number density is increased. 
We do not plot the correlation function for halos below 
rp « 50/i“^ kpc since the halo finder we use does not find 
halos closer than this separation. 

Gonnecting with previous work, we show in the right 
panel of Figure |S1 the halo occupation, {N{M)), for the 
halo samples which correspond to the z ~ 4 Subaru 
samples. Increasingly brighter LBG samples reside in 
preferentially larger mass distinct halos. Furthermore, 
since {N{M)) < 2 for M < our L—14iax 

model implies that the majority of these highest-redshift 
LBGs live alone in distinct dark matter halos. This pic¬ 
ture is in qualitative agreement with previous analysis 
of the LBG clustering in the fram ework of halo model 
llLee et al.ll2005t lOuchi et alJl2005r . 

Figure El compares the clustering of LBGs from 
GOODS data at z ~ 4 (top panels) and z ~ 5 (bot¬ 
tom panels) to the clustering of halos at similar redshifts 
and number densities. Here again the agreement is very 
good, even at z ^ 5. The model somewhat overpredicts 
the clustering of the fainter sample at z = 4. In this plot 
we also show the effect of a contaminated LBG sample, 
i.e. we assume that 10% of identified LBGs are actually 
interlopers with a random distribution on the sky. The 
effect of contamination scales uj{6) by (1 — c)^, where c 
is the contamination fraction. A 10% contamination re¬ 
duces the angular correlation function by « 20% and re¬ 
sults in better agreement for the fainter GOODS samples 
(right panels) which plausibly have higher contamination 
than the brighter samples. We note that 10% contami¬ 
nation is likely an upper limit and was included here for 
illustration purposes. The detailed spectroscopic follow¬ 
up of these high redshift samples have yet to yield direct 
accurate estimates of contamination fractions. 

The straightforward implication of the presented com¬ 
parisons is that the data is consistent with, and one 
could argue supports, the picture in which most LBGs 
are the central galaxies in their host halos with lumi¬ 
nosity tightly related to the halo circular velocity and 
mass. Most LBGs have no neighbors within the same 
halo. However, a fraction of them do and it is this frac¬ 
tion that is responsible for the strong upturn in the corre- 


TABLE 1 

Summary of Samples 


Data 

Defn.“ 

Sb 

n/10-^ 

Mpc“® 

V ^ 

y max 

km 

SDSS 

Mr < -18 

0.04 

27.0 

110 

SDSS 

Mr < -19 

0.06 

15.0 

130 

SDSS 

Mr < -20 

0.06 

6.0 

180 

SDSS 

Mr < -21 

0.15 

1.1 

310 

DEEP2 

Mb < -19.0 

0.87 

13.0 

150 

DEEP2 

Mb < -19.5 

0.92 

8.4 

170 

DEEP2 

Mb < -20.0 

0.98 

4.9 

200 

DEEP2 

Mb < -20.5 

0.99 

2.5 

250 

Adelberger 

UnGR colors 

2.9 

4.0 

207 

Subaru 

i < 25.0 

4 

0.8 

265 

Subaru 

i < 25.5 

4 

1.9 

225 

Subaru 

i < 26.0 

4 

3.8 

191 

Subaru 

i < 26.5 

4 

6.4 

174 

GOODS 

mz < 26.0, S-drop 

4 

2.7 

205 

GOODS 

< 26.5, R-drop 

4 

4.5 

185 

GOODS 

rriz < 26.0, U-drop 

5 

1.5 

200 

GOODS 

mz < 26.5, U-drop 

5 

2.6 

180 


“Absolute SDSS and DEEP2 magnitudes are in units of M — 
51og(h). 

'^Mean redshift of sample; for Subaru and GOODS data, z is 
based on Monte Carlo simulations of artificial LBGs. 

“Imaxsuch that n(> Kmax) = Uganipie for halos at the simula¬ 
tion output closest to z. 


lation function at small scales. By accurately reproduc¬ 
ing both the small-scale upturn in uj{d) and the large- 
scale clustering, our model accurately predicts not only 
the correct distinct halos to associate with LBGs (the ‘2 
halo term’ in halo model jargon) but also the number of 
LBGs within a distinct halo (the corresponding ‘1 halo 
term’). 

5. DEPENDENCE OF HALO OCCUPATION ON NUMBER 
DENSITY AND REDSHIFT 

In this section we explore the redshift and number 
density dependence of the halo occupation distribution 
(HOD), the key ingredient of the halo model, in our 
Fmax —L assignment scheme. We compare our results 
to previous studies which fit the halo occupation to ob¬ 
servational clustering data. The HOD, or the probabil¬ 
ity distribution for a distinct halo of mass M to host 
N galaxies, provides a relatively clean and simple frame¬ 
work for interpreting clustering data. The two-point cor¬ 
relation function in the halo model depends on the first 
and second moments of the HOD {{N) and {N{N — 1)), 
respectively). 

The first moment o f the HOD can be sep arated into 
two components le.g.. iKravtsov et 'mi2004bll : 

{N{M)) = (iV(M))een + (iV(M)),at, (5) 

where {N{M))cen and {N{M))sa.t are the number of cen¬ 
tral and satellite galaxies (i.e., subhalos), respectively. 
The central galaxy term is a step function rising to 
(N(M) )rpn = 1 above a minim um distinct halo mass, 
Mmin. iKravtsov et alJ ll2004bD show that distribution 
of the subhalo occupation at fixed M is well-modeled 
by a Poisson distribution; this continues to hold well in 
the simulations presented here. Although some authors 
model the cutoff at Mmin to take into account scatter be¬ 
tween the observable and host halo mass, we model this 
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term as a step-function for simplicity, and note that this 
simplification will only impact the Mmin parameter. 

The satellite term, {N{M))sat, can be described by a 
power-law fun ction, (N(M))^„t (x M°‘ for l arge distinct 
halo m asses. iKravtsov et al.l l|2004bll and iZheng et al.1 
lj200!T|l find that a « 1 for subhalos and galaxies identi¬ 
fied in cosmological simulations within massive distinct 
halos, over a wide range of number density thresholds 
(corresponding to luminosity threshold) cuts. For small 
distinct halo masses. (N(M))s^t “rol ls-off” faster than 
the power law. IKravtsov et alJl|2004bfl describe the “roll¬ 
off” by: 

(lV(M)),at = (M/Mi-C)“. (6) 

iTinker et al.l ||200 !tI 1 propose an alternative exponential 
form for the “roll-off” : 

(iV(M)),,t = ^exp(-^), (7) 

where Mcut and Mi are free parameters, this time with 
the asymptotic slope fixed at a = 1. In fact, we find that 
Mcut and Mi are tightly correlated such that 

log(A/cut) = 0.761og(Mi) + 2.3, (8) 

for the full range of number densities and redshifts ex¬ 
plored here, thereby reducing Equation [T] to a one- 
parameter function. As shown in Figure I1 1 1 this form 
is an excellent fit to the number of satellite halos in our 
simulations over a large range of redshifts and number 
densities. We have explicitly checked the asymptotic 
slopes of {N{M))sat where possible and find that they 
are consistent with a = 1, in agreement with previous 
theoretical predictions (IKravtsov et alJl2004bll . 

The exponential functional form of Equation 13 is 
preferable to Equation El because for a given choice of 
parameters Mi and C, {N{M))sat can become negative 
using the latter. For example, best-fit values for these 
parameters for the z ~ 0, n = 20.0 x 10“^ Mpc“^ 
sample (triangles in the top left panel of Figure ITTll im¬ 
plies that (iV(M))sat will be negative for M ~ 2 x 10^^ 
h~^MQ, which is not in agreement with the results of 
simulations. For this reason we choose to use Equation 
[3 to describe the functional form of (7V(M))sat- 

Note that the presence of the “roll-off” is important, 
and neglecting it ca n bias the best-fit HOD parameters. 
iZehavi et al.l il200,^ use 

(iV(M)),at = (M/Mi)“ (9) 

to model the clustering of SDSS galaxies at z ~ 0 and 
find that a increases for galaxies of increasing luminos¬ 
ity. However, in the presence of the roll-off at small M, 
brighter samples will have on average fewer satellites and 
hence {N{M))sat will be dominated by the steeper ex¬ 
ponential cut-off. A power-law fit is hence expected to 
yield artificially steep slopes. Indeed, when fitting a pure 
power law without this roll-off to (7V(M))sat in the sim¬ 
ulations, we find slopes closer to a « 1.5. 

In Table|3we list the best-fit HOD parameters for halos 
in our simulations over a wide range of number densities 
and redshifts. The corresponding halo occupations are 
plotted in Figure cni while Figure El shows the contri¬ 
bution from subhalos (satellites) alone. 

With the HOD in hand we can derive several useful 
quantities. The average number of halos per distinct 


TABLE 2 

HOD PARAMETERS AND DERIVED QUANTITIES 


n/lO-3 
h® Mpc“® 

2: 

log(Mniin) 

log(ALi) 

logfAfcut) 

(TV) 

log((TU)) 

1.0 

0 

12.8 

14.0 

13.1 

1.15 

13.6 

1.0 

1 

12.7 

13.7 

12.7 

1.16 

13.2 

1.0 

3 

12.1 

12.0 

12.0 

1.15 

12.5 

1.0 

4 

11.8 

12.4 

12.1 

1.16 

12.2 


4.0 

0 

12.3 

13.5 

12.4 

1.22 

13.4 

4.0 

1 

12.1 

13.3 

12.4 

1.14 

12.9 

4.0 

3 

11.7 

12.8 

12.1 

1.08 

12.2 

4.0 

4 

11.5 

12.5 

11.8 

1.08 

11.8 


8.0 

0 

12.0 

13.1 

12.4 

1.27 

13.4 

8.0 

1 

11.9 

13.0 

12.3 

1.17 

12.8 

8.0 

3 

11.5 

12.6 

11.9 

1.08 

12.0 

8.0 

4 

11.2 

12.3 

11.6 

1.08 

11.7 


20.0 

0 

11.5 

12.7 

12.0 

1.32 

13.4 

20.0 

1 

11.5 

12.7 

11.9 

1.20 

12.7 

20.0 

3 

11.1 

12.4 

11.6 

1.09 

11.8 

20.0 

4 

10.9 

12.1 

11.5 

1.07 

11.5 


Note. — All quoted masses are in units of h ^Mq. 
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Fig. 10.— Evolution of the halo occupation, {N{M)) from z ~ 4 
to 0 at fixed number density (clockwise from top left: n = 1.0, 4.0, 
20.0, and 8.0, in units of Mpc“®). 


halo, 


(TV) 


(jV(M))nhaio(M)dM 

«halo(M)dM 


( 10 ) 


and the average mass of distinct halos (member number 
weighted), 


(M) 


(TV(M))nhaio(M)MdM 

”‘”/M..-halo(M)dM 


( 11 ) 


are useful when considering a ‘typical’ halo in a given 
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sample. These two quantities are included in Table [21 
Note that when computing these quantities we use 
the halo occupation measured directly from simulations 
(shown in Figure [T(1ll rather than the best-ht halo occu¬ 
pation implied by equations Oand 0 

The redshift evolution of {N{M)) for samples of fixed 
number density reveals several interesting trends. The 
minimum mass of the sample, or the location of the step 
in the mean occupation number, is decreasing with in¬ 
creasing redshift, reflecting evolution of the halo mass 
function. At higher redshifts, the halos hosting multi¬ 
ple satellites become more rare, so that in most samples 
the high-M power-law tail of {N{M))sat is not present. 
This trend is apparent in Figure lTTI which shows the red¬ 
shift evolution of {N{M))sat for halo samples of differ¬ 
ent number densities. At the same time, the “shoul¬ 
der” (or the region between the step and the power- 
law tail) in {N{M)) becomes shorter and not as flat at 
higher z (Fig. 11 OH . reflecting the increasing fraction of 
relatively small-mass halos hosting more than a single 
central galaxy. We choose to characterize the extent of 
the shoulder via the ratio Mi/Mjain- As can be seen by 
consulting Table (21 this ratio (and hence the shoulder) 
decreases both with decreasing number density and in¬ 
creasing redshift. As we will discuss more thoroughly 
below, we believe that it is the extent of the shoulder 
which is primarily responsible for the upturn in the cor¬ 
relation function on small scales. Hence samples with 
a lower number density and/or at higher redshift (and 
hence a smaller shoulder) should have a more significant 
upturn in ojpirp) on small scales, as is observed. 

The HOD parameters we derive are in go od qualitative 
agreem ent with observations. At z ~ 0, iZehavi et al.l 
IMS) fit diffcrGiit functionSjl form for 
(namely Eqn. o, and hence a direct comparison is com¬ 
plicated. However, Mmin is unaffected by the functional 
form for {N{M))sat, and here the values quoted in Ta - 
ble|21are in excellent agreement with iZehavi et all ll2005ll . 
The trend of Mi to decrease with increasing number den¬ 
sity is also in agreement with the derived results from 
observations. 

A direct comparison to the HOD derived in the lit¬ 
erature for LBGs is again complicated. The parame¬ 
ters given in Table (21 are predicted by our simulations, 
while most analyses assume a pure power-law form for 
the overall halo occupation and fit fo r the parame- 

l|Hama,na et al.l 

. The assumed 
power-law form does not describe HOD in our simula¬ 
tions well. This is an example illustrating that HOD 
fitting is in general not unique — different functional 
forms can simultaneously match the number density 
and 2-point clustering of a sample. Because derived 
quantities such as {N) and (M) are entirely dependent 
on the form of the HOD, it is critical to realize that 
these derived quantities are also not unique. Despite 
the different forms for {N{M)), we find here similar 
qualitative trends for the average distinct halo mass of 
LBGs, namely that fainter samples (corresponding to 
higher number den sitie s) have smaller ( M), as is seen in 
lOuchi et alJ ll2005ll and^g^^yj|2Q]2^jb As emphasized 
bv iBullock. Wechsler. fc Somervillel l|2002ll . this degener¬ 
acy may be partially broken by looking at the smallest 
scale clustering information available (typically parame¬ 


ters in the framework of the halo model 




10 “ 10 ‘" 10 “ 10 “ 10 “ 10 “ 10 “ 10 “ 10 “ 
Mass (h'‘ M@) 


Fig. 11.— {N{M))sat as a function of number density and red¬ 
shift. Diamonds designate halo samples at n = 4.0 x 10“^ 
Mpc~^, triangles are for samples with n = 20.0 x 10~^ Mpc~^. 
The dashed lines corresponds to Eqn. [3 with the best-fit parame¬ 
ters listed in Table [3 


terized as the close-pair counts); a comparison with the 
most recent data of this sort from z ^ 0 — 1 is being 
pursued for a similar model to ours by Berrier et al (in 
preparation). 

As was shown in m selecting subhalos according 
to results in an enhanced clustering signal com¬ 
pared to using j especially on small scales and at 

low redshifts, because there are more subhalos in the 
^max selected sample. To illustrate this directly, we fit 
Eqn. 0 to satellites selecting in these two different ways. 
We find that the effect of using ^max i inso¬ 

far as there are more satellites at fixed distinct halo 
mass, is small for low masses, and gradually increases 
toward higher masses. In other words, the difference in 
{N{M)) is small near the exponential cutoff while 
selected samples show steeper {N{M)) with increasing 
M for higher masses. 

To understand this trend consider the following. In 
order for there to exist a sizable difference between 
and T/Jax > n subhalo must orbit within a distinct 
halo for a sufficient amount of time such that tidal strip¬ 
ping can reduce Knax ■ Dynamical friction, which de¬ 
pends on the relative masses of subhalos to distinct ha¬ 
los, will be much stronger in lower distinct halo masses 
because there exists a minimum halo mass corresponding 
to our Knax threshold. In other words, as one considers 
smaller distinct halo masses for a given Knax threshold, 
the typical subhalo mass cannot be arbitrarily close to 
the host mass because such halos would suffer efficient 
dynamical friction and would merge quickly without 
significant evolution of Knax • At higher distinct halo 
masses, dynamical friction is not efficient for the typical 
subhalo since now the typical subhalo is much less mas¬ 
sive than the distinct halo, and hence a subhalo can exist 
within a distinct halo for a sufficient amount of time for 
tidal stripping to have a significant effect on Knax ■ 
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6. DISCUSSION 

Following the reliable identification of subhalos in N- 
body simulations, there has been a concerted effort to 
identify such subhalos with observed galaxies. A per¬ 
sistent problem with this approach is that subhalos se¬ 
lected using the present-day mass or circular velocity 
tend to have radial distributions within their larger dis¬ 
tinct halos that are shallower than those exhibited by 


observed galaxies 

Gao et al.ll2004 IDiemand et alJl2004 

Nagai & Kravtsov 

1200,'ll iMaccio et alJ l200,'jll . However, 

Nagai & Kravtsov 

1200^ find that the observed radial 


profile of galaxies can be reproduced by selecting subha¬ 
los using stellar mass, rather than the present day total 
mass. We argue that the best proxy for the stellar mass 
in dissipationless simulations is the current maximum cir¬ 
cular velocity for distinct halos and the circular velocity 
at the time of accretion, U^ax ; for subhalos. 

We have shown that this simple, non-parametric model 
which relates galaxy luminosities to the maximum cir¬ 
cular velocity of dark matter halos, Umax, by preserv¬ 
ing the observed galaxy luminosity function, accurately 
reproduces the observed luminosity-dependent cluster¬ 
ing of galaxies over projected separations O.I < rp/{h~^ 
Mpc) < 10.0, from z ~ 5 to 0 ^ 0. The only assumption 
of the model is that there is a monotonic relation between 
galaxy luminosity and halo circular velocity. To assign 
luminosities we use the maximum circular velocity at the 
time of accretion, U,^™ j for subhalos, and the maximum 
circular velocity at the time of observation, Ujj5°x , for 
distinct halos. This ingredient is crucial for accurately 
reproducing the observed clustering of galaxies at z ~ 0. 

The success of the model has several important im¬ 
plications. First, our results indicate that dissipation¬ 
less simulations do not suffer from significant overmerg¬ 
ing because they reproduce both the amplitude and 
shape of the correlation function and the underlying 
HOD quite accurately. If our dissipationless simulations 
were missing a significant fraction (> 40%) of objects in 
groups and clusters, as suggested by so me recent studies 
llGao et ^12004 iDiemand et al.lU004ll . both the ampli¬ 
tude and shape of the correlation function in these sim¬ 
ulations would be grossly incorrect. We also note that 
recent analysis of the observed galaxy-mass cross corre¬ 
lation function disfavors the existence of a large popu- 
lation of satellites withou t associated dark matter halos 
ijMandelbaum et alJl200^. 

In order to assess the quantitative impact of a fraction 
of “orphan” gala xies — satellite galaxies with no iden¬ 
tifiable subhalo llGao et aLll2nn4|l — in our model, we 
perform the following test. We artificially increase the 
fraction of subhalos in such a way that the subhalo occu¬ 
pation function, {N(M))sat increases in amplitude while 
maintaining the same shape. In order to match the ob¬ 
served galaxy luminosity function, we simulatenously in¬ 
crease the overall Umax threshold for a given sample, such 
that the number density does not change. We find that 
an orphan fraction > 10% is inconsistent with the cor¬ 
relation functions of the SDSS galaxy samples at z ^ 0 
for all luminosities (for this cosmological model and the 
assumption of no scatter between L and U^ax )• 

Our results also imply that the central assumption of 
the luminosity assignment model — the tight, mono¬ 
tonic relation between galaxy luminosity and halo cir¬ 


cular velocity — likely exists for the o bserved galaxies . 
Such a relation is indeed expected (e.g., iMo et ai]ll998ll 
for isolated galaxies, but our results indicate that this 
is true globally for galaxies of different types and for 
a wide range of redshifts. We argue that for subhalos 
the dissipation should result in a centrally condensed, 
tightly-bound stellar system which would stabilize the 
halo circular velocity against tidal heating. Stellar mass 
and luminosity of galaxies should therefore correlate with 
the circular velocity of subhalos at the time they are 
accreted, before significant tidal evolution takes place. 
Given that we match luminosity at a particular epoch 
to the circular velocity at different epochs (the epoch of 
observation for distinct halos, and the epoch of accretion 
for subhalos), a subtle implication of our scheme is that 
the relation between luminosity and Umax does not evolve 
strongly with time — a result which may have been an¬ 
ticipated by the lack of scatter in the Tully-Fisher rela¬ 
tionship in different environments. 

The corollary is then that the clustering of a particular 
galaxy sample is largely determined by the clustering of 
halos and subhalos that host the galaxies. Clustering of 
halos and subhalos is governed by gravitational dyna mics 
fe.g.. iKravtsov fc KlvDinlll999l: IZentner et ani200fi[ and 
discussion therein), while the particular subset of halos 
that host galaxies in a given sample is determined by 
the relations between observable galaxy properties and 
properties of their host dark matter halos and selection 
criteria used to define the sample. In the case of the 
galaxy luminosity, the relation with halo circular velocity 
appears to be particularly tight. 

The model agreement with the clustering properties 
of the LEG population at 3 < z < 5 is perhaps 
more surprising, given that the LEG selection criteria 
are significantly more complicated than those of the 
luminosity-selected samples at z < I. The success of 
our L—Umax model indicates that the rest-frame LEG 
luminosity is likely to be tightly related to the halo cir¬ 
cular velocity (and hence total mass). This would, in 
turn, suggest that LBGs are fair tracers of the overall 
halo population rather than a special subset of halos, 
such as halos undergoing mergers or collisions. Although 
we have not explicitly investigated the latter scenarios 
(i.e., mergers or collisions), it seems unlikely that the ob¬ 
served LEG clustering would be well matched by a model 
in which the effective duty cycle for the LEGs is small 
and most LEGs are associated with minor mergers or 
small-mass collisions. We note, however, that differences 
in the clustering properties and inferred host masses for 
the so-called “massive halo” model and the “collisional 
starburst” model may have been exaggerated — in the 
former case because the massive subhalos were ignored, 
and in the latter because the efficiency of star formation 
in me rgers was probably overestimated (e.g., iGox et all 

unns). 

A key success of the simulations and our model is the 
correct description of the small-scale upturn in the cor¬ 
relation functio n. At z ~ 0, the up turn is detected for 
bright galaxies llZehavi et al.ll2(i(M) . while the correla¬ 
tion function of faint o bjects is very close to a power 
law llZehavi et aLll^DlTfl . At z > I, the upturn is more 
pronounced and is now unambiguously detected both in 
the DEEP2 data and in LEG samples at z ~ 3 — 5, as 
predicted from the halo model arguments llZhentHl2nn4^ 
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and cosmological simulations llKravtsov et al.l l2004bll . 
The trend towards more pronounced deviation from the 
power-law form of CF at higher z appears to be due to 
a couple different factors. 

The accretion rate of subhalos is substantially higher 
at high z, so that the abundance of subhalos with masses 
close to the threshold mass of the sample, Mmin, is also 
larger on average (see § 0- This makes the “shoul¬ 
der” of the HOD of high-z objects shorter and not as 
flat, increasing the number of small-separation pairs be¬ 
tween central galaxies and satellites and between satel¬ 
lites themselves in small-mass halos. This also makes the 
width of the HOD at masses close to Mmin wider and 
closer to the Poisson distribution. Although the overall 
fraction of subhalos in a given sample actually decreases 
with increasing redshift (see Fig. the average mass of 
the distinct halos hosting most of subhalos is decreasing 
as well. The contribution of the pairs of objects within 
the same halo to the 1-halo term of the correlation func¬ 
tion comes then predominantly from the smallest halos in 
the sample and the 1-halo term is considerably “peakier” 
than at lower redshifts, where halos of a wide range of 
masses contribute. We argue therefore that the more 
pronounced upturn in the correlation function at higher 
redshifts and for rarer objects at a given epoch is primar¬ 
ily due to the increased fraction of subhalos in host halos 
with masses closer to the threshold mass of the sample. 

In the present work, we have not included any scat¬ 
ter in matching luminosity and velocity relations, and 
it is interesting that making this approximation works 
so well. Some scatter in the A—Vmax relation should 
certainly be expected given the scatter in the observed 
Tully-Fisher and Faber-Jackson relations. This scatter 
likely depends on band, galaxy properties, and the ex¬ 
tent to which internal dust extinction has been corrected 
for. For our purposes, the main effect of adding scatter 
in this relationship to our model would be at the high¬ 
est luminosities. The insensitivity at low luminosities 
arises because 1) the luminosity function of galaxies is 
rather flat at A < A* so that the scatter results in al¬ 
most equal numbers of objects scattering in and out of 
the sample, and 2) the bias of halos corresponding to 
these luminosity cuts (A < A*), is only a weak func¬ 
tion of mass and luminosity both in o ur simulations and 
in observations llTegmark et alJl2004fl . so changes in the 
mass distribution of galaxies do not affect clustering sig¬ 
nificantly. In contrast, at higher luminosities, both the 
luminosity function and the bias dependence on mass are 
steep, and adding scatter serves to decrease the average 
clusterin g and makes the l umin osity a steeper function 
of mass. iTasitsiomi et al ] illiol find that some scatter 
is required to match the galaxy-mass correlations, but 
the conclusion is reached using IVJax subhalos, which 
we expect to exhibit more scatter with luminosity than 
Knax • There is, however, an interesting tension between 
constraints provided by the observed galaxy-galaxy and 
galaxy-mass correlation functions for the most massive 
halos (cf. Fig. |3J): the former constrains the scatter for 
bright galaxies to be small, while the latter appears to 
require significant amount of scatter. We postpone a full 
discussion of these issues including combined constraints 
to future work. 


7. SUMMARY 

Our main results and conclusions are as follows: 

1. We show that a simple, non-parametric model 
which monotonically relates galaxy luminosities to 
the maximum circular velocity of dark matter ha¬ 
los, Vmax , by preserving the observed galaxy lumi¬ 
nosity function, accurately reproduces the observed 
luminosity-dependent clustering of galaxies over 
projected separations 0.1 < rp/{h~^ Mpc) < 10.0, 
through most of the evolution of the universe from 
z ^ 5 to the present. 

2. The key to the success of the model is our luminos¬ 
ity assignment scheme, in which we use the max¬ 
imum circular velocity at the time of accretion, 

j for subhalos and the maximum circular ve¬ 
locity at the time of observation, j for distinct 
halos. We argue that for subhalos in dissipationless 
simulations, K^ax reflects the luminosity and stel¬ 
lar mass of the associated galaxies better than the 
circular velocity at the epoch of observation, k(((ax 

3. Our simulations, coupled with the above lumi¬ 
nosity assignment scheme, correctly reproduce the 
small-scale deviation of the correlation function 
from the power-law form (the upturn) in observed 
luminosity-selected samples at different redshifts. 
The deviation is, in general, more pronounced at 
high redshifts. We attribute this trend to the in¬ 
creased fraction of satellite galaxies in host halos 
with masses close to the threshold mass of the sam¬ 
ple. 

4. Our luminosity assignment model applied to the 
simulations at z ~ 3 — 5 reproduces the observed 
shape and amplitude of the two-point correlation 
function of Lyman-break galaxies on both large and 
small scales. This suggests that LBGs are fair trac¬ 
ers of the overall halo population at these epochs 
and that their luminosity can be simply related to 
the circular velocity (and hence mass) of their dark 
matter halos. Although our model implies that 
most (~ 90%) of the LBGs are central galaxies in 
distinct halos, a fraction of the LBGs correspond 
to subhalos (i.e., satellites). This fraction is suffi¬ 
ciently high for a significant deviation of the LBG 
correlation function from the power law at small 
scales. 

It is quite remarkable that, after vast improvements in 
observational data over the past decade, the simple pic¬ 
ture for the for mation of galaxies wi t hin virialized dark 
matter halos ilWhite fc ReesI 119781 iFall fc EfstathioiJ 
IT^ . whereby galaxy properties such as luminosity are 
tightly coupled to the properties of their host halos, pro¬ 
vides a good description of the observed clustering trends 
of galaxies over nearly the entire age of the Universe, on 
both large and small scales. 

As we discussed above, our A—14iax model has nat¬ 
ural implications for the Tully-Fisher relation and dy¬ 
namical mass-to-light ratio, and the evolution in these 
quantities. A fruitful followup to the study presented 
here would be to compare the predictions of this model 
to local and to (forthcoming) DEEP2 measurements of 
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the relation. The comparison of model predictions to ob¬ 
served galaxy-galaxy and galaxy-mass clustering can also 
constrain the amount of scatter in the L—Umax relation 
for bright galaxies, and readily predicts the conditional 
luminosity function of galaxies (luminosity function as a 
function of distinct halo mass), and its evolution. Our lu¬ 
minosity assignment scheme can also be extended by sep¬ 
arately modeling central and satellite galaxies, which ap¬ 
pear to be two distinct po pulations with different origins 
and observable properties iVale fc Ostrikeil2004HGo oravl 
l2005alhl : iWeinmann et alJ|2005HHansen et al.ll2005ll . 

The simulations were run on the Columbia machine 
at NASA Ames and on the Seaborg machine at NERSC 
(Project PI: Joel Primack). First and foremost, we would 
like to thank Anatoly Klypin for running these simu¬ 
lations and making them available to us. We are also 
indebted to Brandon Allgood for providing the merger 
trees that were used to measure U^ax • We thank Kurt 
Adelberger, Alison Coil, Mauro Giavalisco, Kyoungsoo 
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